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ABSTRACT 

We present a candidate for the most distant galaxy known to date with a photometric redshift 
z = 10.7^0 4 (95% confidence hmits; with z < 9.5 galaxies of known types ruled out at 7.2-<j). This 
J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing 
and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due 
to strong gravitational lensing by the galaxy cluster MACS J0647. 7+7015 at z = 0.591. The images 
are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB 
(-0.15/iJy) in the WFC3/IR F160W filter ( -1. 4-1. 7/im) where they are detected at >12-cr. Ah three 
images are also confidently detected at ^6-<j in F140W (~1.2-1.6/im), dropping out of detection from 
15 lower wavelength HST filters (~0.2-1.4/im), and lacking bright detections in Spitzer/IRAC 3.6/im 
and 4.5/im imaging (~3.2-5.0/im). We rule out a broad range of possible lower redshift interlopers, 
including some previously published as high redshift candidates. Our high redshift conclusion is 
more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH 
observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and 
MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower 
redshift s. This would suggest that low luminosity galaxies could have reionized the universe. However 
given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop 
off in number counts at z > 10 suggested by field searches. 

Subject headings: early universe - galaxies: high redshift - galaxies: distances and redshifts - 
galaxies: evolution - gravitational lensing: strong - galaxies: clusters: individ- 
ual (MACSJ0647.7+7015) 
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1. INTRODUCTION 

Current models of structure formation suggest that the 
first galaxies formed at z > 10 when the universe was 
< 470 million years old (Wise & Abel 2007; Wise et al. 
2008; Greif et al. 2008, 2010; and for recent reviews, see 
Bromm & Yoshida 2011 and Dunlop 2012). Observations 
may be closing in on these first galaxies with one z ~ 
10 candidate detected in the Ultra Deep Field (UDFj- 
39546284; Bouwens et al. 2011a) and another strongly 
lensed by a galaxy cluster (MACS 1 149- JDl; Zheng et al. 
2012). 

Intriguingly, the number density of z ~ 10 galaxies 
detected in unlensed fields is several times lower than 
predicted based on extrapolations from lower redshifts, 
assuming a luminosity function with one or more param- 
eters evolving linearly with redshift (Bouwens et al. 2008, 
2011a; Oesch et al. 2012a). This suggests that the star 
formation rate density built up more rapidly from z ~ 10 
to 8 than it did later between z ~ 8 and 2. This is con- 
sistent with some theoretical predictions (Trent i et al. 
2010; Lacey et al. 2011). However, Robertson & Ellis 
(2012) suggest such a sharp drop off would be in tension 
with z < 4 GRB rates as correlated with star formation 
rate density and extrapolated to higher redshifts. Di- 
rect detections and confirmations of z > 10 galaxies are 
required to more precisely constrain the star formation 
rate density at that epoch. 

The observed luminosity functions at z ~ 7 and 8 fea- 
ture steep faint end slopes of a ~ — 2 (Bouwens et al. 
2011b; Bradley et al. 2012b), steeper than at lower red- 
shifts, a trend consistent with model predictions (Trenti 
et al. 2010; Jaacks et al. 2012). If these luminosity func- 
tions can be extrapolated to 2: > 10, then low lumi- 
nosity galaxies (Mjjy fainter than —16 AB) could have 
reionized the universe (Bouwens et al. 2012a; Kuhlen & 
Faucher-Giguere 2012), assuming a sufficient fraction of 
their UV photons could escape their host galaxies to the 
surrounding medium (see also Conroy & Kratter 2012). 
Otherwise, a more exotic source of reionizing energy may 
have been required, such as self-annihilating dark matter 
(locco 2010; Natarajan 2012). 

Reionization was likely well underway by 2; > 10 but 
with over half the universe still neutral (Robertson et al. 
2010; Pandolfi et al. 2011; Mitra et al. 2012). Improving 
our understanding of the early universe and this phase 
change is one of the pressing goals of modern cosmology. 

Observations with the Wide Field Camera 3 (WFC3; 
Kimble et al. 2008) installed on the Hubble Space Tele- 
scope (HST) have significantly advanced our understand- 
ing of the z^l universe, over 13 billion years in the past. 
The Ultra Deep Field and surrounding deep fields have 
yielded over 100 robust z >1 candidates as faint as 29th 
magnitude AB (Bunker et al. 2010; Labbe et al. 2010; 
Bouwens et al. 2011b; Oesch et al. 2012a). Analyses of 
wider space-based surveys such as CANDELS (Grogin 
et al. 2011; Koekemoer et al. 2011), BoRG (Trenti et al. 
2012), and HIPPIES (Yan et al. 2011) have helped fill 
out the brighter end of the luminosity function (Oesch 
et al. 2012b; Bradley et al. 2012b; Yan et al. 2011). 

Of the handful of z > 7 galaxies spectroscopically con- 
firmed to date, most have been discovered in even wider 
near-infrared surveys carried out with the ground-based 
telescopes Subaru (Shibuya et al. 2012; Ono et al. 2012), 



VLT (Vanzeha et al. 2011), and UKIRT (Mortlock et al. 

2011) . Surveys with the VISTA telescope are also be- 
ginning to yield high redshift candidates (Bowler et al. 

2012) . 

Complementary to these searches of "blank" fields are 
searches behind strongly lensing galaxy clusters (Kneib 
et al. 2004; Bradley et al. 2008; Bouwens et al. 2009; 
Bradac et al. 2009; Maizy et al. 2010; Richard et al. 2011; 
Hah et al. 2012; Bradley et al. 2012a; Zitrin et al. 2012a; 
Wong et al. 2012; Zackrisson et al. 2012; Bradac et al. 
2012; and for a recent review, see Kneib & Natarajan 
2011). The drawbacks of lensed searches are reduced 
search area in the magnified source planes and some un- 
certainty in the estimate of that search area introduced 
by the lens modeling. But the rewards are galaxies which 
are strongly magnified, often by factors of 10 or more. 
Lensed searches are significantly more efficient in yield- 
ing high-redshift candidates bright enough for spectro- 
scopic confirmation, including A1703-zD6 (Bradley et al. 
2012a) at z = 7.045 (Schenker et al. 2012). 

The Cluster Lensing and Supernova survey with Hub- 
ble (CLASH; Postman et al. 2012) is a large Hubble pro- 
gram imaging 25 galaxy clusters in 16 filters including 
five in the near-infrared (0.9-1.7/im). Five of these, in- 
cluding MACSJ0647.7+7015 {z = 0.591, Ebeling et al. 
2007), were selected on the basis of their especially strong 
gravitational lensing power as observed in previous imag- 
ing, with the primary goal of discovering highly magni- 
fied galaxies at high redshift. To date, some of the more 
notable strongly lensed galaxies found in CLASH include 
a doubly-imaged galaxy with a spectroscopic redshift 
z = 6.027 observed at mag ~ 24.6 (Richard et al. 2011) 
which is possibly ~800 Myr old (although see Pirzkal 
et al. 2012); a quadruply-imaged z ~ 6.2 galaxy ob- 
served at 24th magnitude (Zitrin et al. 2012a); and the 
z ~ 9.6 candidate galaxy MAC SI 149- JDl observed at 
mag - 25.7 (Zheng et al. 2012). The z - 9.6 candidate is 
strongly lensed by MACSJ1149.6-I-2223, another CLASH 
cluster selected for its high magnification strength. 

Here we report the discovery of MACS0647-JD, a can- 
didate for the earliest galaxy yet detected at a red- 
shift z = 10.7^0.4 (95% confidence), just 427^2? mil- 
lion years after the Big Bang. It is strongly lensed by 
MACSJ0647.7+7015, yielding three multiple images ob- 
served at F160W AB mag -25.9, 26.1, and 27.3, magni- 
fied by factors of —8, 7, and 2. The brightest image is 
similar in fiux to MACS1149-JD1 (F160W mag - 25.7) 
at 2: = 9.6 ±0.2 (68% confidence) and roughly 15 times (3 
magnitudes) brighter than the z = 10.3 ± 0.8 (68% con- 
fidence) candidate in the UDF (Bouwens et al. 2011a). 

MACS0647-JD is a J-dropout as all three lensed im- 
ages are securely detected in F160W and F140W but 
drop out of detection in the J-band F125W and all 14 
bluer HST filters. We show this photometry is most likely 
due to the Lyman-a break redshifted to — 1.46/im at 
z ^ 11. This Lyman dropout technique (Meier 1976; Gi- 
avalisco 2002) pioneered by Steidel et al. (1996) at z — 3 
has been used with a high success rate to identify high- 
redshift candidates later spectroscopically confirmed out 
to z — 7. However care must be taken not to confuse 
dropouts with intrinsically red (evolved and/or dusty) 
galaxies at intermediate redshift (Schaerer et al. 2007; 
Dunlop et al. 2007; Chary et al. 2007; Capak et al. 2011; 
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Table 1 

Observed Filters and Integration Times 



Filter wavelength^ exposure 
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F475W 


0.47/im 


2248 


sec 


F555W 


0.54/im 


7740 


sec 


F606W 


0.59/im 


2064 


sec 


F625W 


0.63/im 


2131 


sec 


F775W 


0.77/im 


2162 


sec 


F814W 


0.81/im 


12760 


sec 


F850LP 


0.90/im 


4325 


sec 


F105W 


1.06/im 


2914 


sec 


FllOW^ 


1.15/im 


1606 


sec 


F125W 


1.25/im 


2614 


sec 


F140W 


1.39/im 


2411 


sec 


F160W 


1.54/im 


5229 


sec 


IRAC chl 


3.55/im 


18000 


sec 


IRAC cli2 


4.50/im 


18000 


sec 



^ Effective "pivot" wavelength (Tokunaga 
& Vacca 2005). 

2 Visit A2 only, excluding visit A9 (§3.1). 



Boone et al. 2011; Hayes et al. 2012). In the case of 
MACS0647-JD, we show it is extremely difficult, if not 
impossible, for low redshift interlopers to reproduce the 
observed colors, especially the J125 — i^ieo ^ 3 magni- 
tude break of MACS0647-JD1. We also test our analysis 
method by reanalyzing previously published J-dropouts 
which later proved to be at intermediate redshift. Our 
Bayesian photometric redshift (BPZ, Benitez 2000; Coe 
et al. 2006) analysis correctly shows intermediate redshift 
solutions are preferred for those objects, while preferring 
a higher redshift for MACS0647-JD. 

We describe our HST and Spitzer observations in §2 
and present photometry in §3. We derive the photomet- 
ric redshift in §4 and consider a wide range of possible 
interlopers in §5. We present our gravitational lensing 
analysis in §6. In §7 we derive physical properties of 
MACS0647-JD based on additional photometric analy- 
sis. In §8 we compare our observed number density of 
z ^ 11 galaxies to that expected, and we constrain the 
z > 9 star formation rate density. Finally, we present 
conclusions in §9. 

Where necessary we assume a concordance ACDM cos- 
mology with h = 0.7, Vtrn = 0-3, ^^a = 0.7, where 
Hq = 100 h km s-^ Mpc-^ In this cosmology 1" = 3.93 
kpc at z = 10.8 and 6.62 kpc at the cluster redshift 
z = 0.591. 

2. OBSERVATIONS 

As part of the CLASH program, HST observed the 
core of MACSJ0647.7+7015 (Fig. 1) during 19 orbits 
spread among eight different visits between October 
5 and November 29, 2011 (General Observer program 
12101). Imaging was obtained with the Wide Field Cam- 
era 3 (WFC3; Kimble et al. 2008) and Advanced Camera 
for Surveys (ACS; Ford et al. 2003) in 15 filters spanning 
0.2-1.7/im including five near-infrared WFC3/IR filters 
spanning 0.9-1.7/im. These datasets were supplemented 
by prior ACS imaging obtained in the F555W (0.56/im) 
and F814W (0.81/im) filters to total depths of - 3.5 and 
5.5 orbits, respectively (GO 9722 P.I. Ebeling; GO 10493, 



10793 P.I. Gal- Yam). These observations are detailed in 
Table 1. 

We processed the images for debias, (super)flats, and 
darks using standard techniques, then co- aligned and 
combined them to a scale of 0.065" /pixel; see Koeke- 
moer et al. (2002, 2011) for further information on the 
astrometric alignment and drizzle algorithms that were 
used and Postman et al. (2012) for specific details on 
their implementation in CLASH. We also produced in- 
verse variance maps (IVMs) based on the observed sky 
level, identified cosmic rays, detector flat field, readnoise, 
dark current, and bad pixels. These IVMs may be used 
to estimate the level of uncertainty in each pixel before 
accounting for correlated noise and any Poisson source 
noise. 

Imaging at longer wavelengths was obtained by the 
Spitzer Space Telescope with the InfraRed Array Cam- 
era (IRAC; Fazio et al. 2004) chl (3.6/im) and ch2 
(4.5/im) with total exposure times of 5 hours at each 
wavelength (Program 60034, P.I. Egami). These obser- 
vations were divided into two epochs separated by ~5.5 
months (Nov. 10, 2009 and Apr. 23, 2010). We com- 
bined the Basic Calibrated Data (BCD) using MOPEX 
(Makovoz & Khan 2005) to produce mosaicked images. 

As of July 2012, CLASH had obtained 16- 
band HST observations for 17 clusters, includ- 
ing MACSJ0647.7+7015 and three other "high- 
magnification" strong lensing clusters, as given in 
Table 3. We searched for high redshift galaxies in 
the WFC3/IR fields of view of all 17 of these clusters 
(Bradley et al. 2012c, in preparation). 

Out of ~20,000 detected sources, we identified 
MACS0647-JD (Fig. 2) as having an exceptionally high 
photometric redshift (§4). Our selection was based on 
SED (spectral energy distribution) fitting as used in some 
previous high-redshift searches (e.g., McLure et al. 2006; 
Dunlop et al. 2007; Finkelstein et al. 2010). We did not 
impose specific magnitude limits, color cuts, or other de- 
tection thresholds on our selection as in other works (e.g.. 
Bunker et al. 2010; Yan et al. 2011; Bouwens et al. 2012b; 
Oesch et al. 2012b). 

3. PHOTOMETRY 
3.1. HST Photometry 
3.1.1. Photometric Analysis 

We used SExtractor version 2.5.0 (Bertin & Arnouts 
1996) to detect objects in a weighted sum of all five HST 
WFC3/IR images. Along the edge of each object, SEx- 
tractor defines an isophotal aperture consisting of pixels 
with values above a detection threshold. We set this 
threshold equal to the RMS measured locally near each 
object. Isophotal fluxes (and magnitudes) are measured 
within these isophotal apertures. SExtractor derives flux 
uncertainties by adding in quadrature the background 
RMS derived from our inverse variance maps and the 
Poisson uncertainty from the object flux. 

Since our images are drizzled to a 0.065" pixel scale, 
which is 2-3 times smaller than the WFC3 PSF, the re- 
sulting images contain significant correlated noise. The 
weight maps produced by drizzle represent the expected 
variance in the absence of correlated noise. To account 
for the correlated noise, one may apply a correction fac- 
tor as in Casertano et al. (2000). 
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Figure 1. Lenstool strong lensing mass model of M ACS J0647. 7+7015 and multiply-imaged galaxies as identified in this work using the 
Zitrin et al. (2009) method, including two strong lensing systems identified in Zitrin et al. (2011a). Each strongly lensed galaxy is labeled 
with a number and color-coded by redshift (scale at bottom right). Letters are assigned to the multiple images of each galaxy. Dashed 
circles indicate predicted locations of counterimages not unambiguously identified. Overplotted are critical curves from our Lenstool model 
indicating thin regions of formally infinite magnification for background galaxies at z = 2.0 (cyan), 3.5 (green) and 11.0 (red). Mirror 
images of galaxies straddle these critical curves. The Hubble color image was produced using Trilogy (Coe et al. 2012) and is composed of 
ACS and WFC3/IR filters as given at top right. 



Previous authors have also noted that SExtractor 
tends to underestimate flux uncertainties (Feldmeier 
et al. 2002; Labbe et al. 2003; Gawiser et al. 2006) by 
as much as a factor of 2-3 (Becker et al. 2007). In this 
work, we obtained empirical measurements of the flux 
uncertainties using the following method which also cap- 
tures the effects of correlated noise. 

SExtractor has the ability to measure the local back- 
ground within a rectangular annulus (default width 24 
pixels) around each object. We constructed a rectangle 
of the same size, but rather than calculate the RMS of 
the individual pixels, we obtained samples of the back- 
ground flux within this region using the isophotal aper- 
ture shifted to new positions. In other words, we moved 
the isophotal aperture to every position within this rect- 
angle, sampling the flux at each position. We discarded 
measurements for which the aperture includes part of 
any object, as we are interested in measuring the back- 
ground flux. Finally we measured the RMS of these mea- 
surements and added in quadrature the object's Poisson 
uncertainty to obtain the total flux uncertainty for that 
object. We found this technique indeed yielded larger 
flux uncertainties than reported by SExtractor, typically 
by factors of 2-3 in the WFC3/IR filters and by lower 
factors in ACS and WFC3/UVIS. 

We also used this method to determine object fluxes. 
The mean of the flux measurements in the nearby aper- 
tures was adopted as the local flux bias, which we sub- 



tracted from the flux measurement in the object itself. 
We found this yielded photometry very similar to that 
obtained using SExtractor, agreeing well within the pho- 
tometric uncertainties. While we used this photometry 
for all subsequent analyses, we also verified that our de- 
rived photometric redshifts did not vary significantly (af- 
ter excluding the FllOW second epoch exposures; §3.1.3) 
if we instead utilized photometry derived directly from 
SExtractor. 

We corrected for Galactic extinction of E{B — V) = 
0.11 in the direction of MACSJ0647.7+7015 as derived 
using the Schlegel et al. (1998) IR dust emission maps. 
For each filter, the magnitudes of extinction per unit 
E{B — V) are given in Postman et al. (2012, their Table 
5). (These values should be ~10% lower in the NUV and 
optical according to Schlafly & Finkbeiner 2011.) This 
extinction reddens the observed colors at the few percent 
level in the near-IR. Thus the effects on the J-dropout are 
negligible. The extinctions range from 0.05 to 0.11 mag 
in the WFC3/IR images; 0.16 to 0.46 mag in the ACS 
images; and 0.50 to 0.83 mag in WFC3/UVIS. We note 
the extinction may be somewhat uncertain due to patchy 
galactic cirrus in the direction of MACSJ0647. 7+7015. 

3.1.2. Photometric Results 

Our resulting 17-band HST photometry is given in Ta- 
ble 2 and Fig. 4. All three J-dropouts are detected at 
>10-cr in F160W, >6-a in F140W, and <3-cr in ah other 
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Table 2 

Coordinates, Observed Filters, and Photometry of the J-dropouts 



JDl 



JD2 



JD3 



JDl + JD2 + JD3^ 



R.A. (J2000) 

Decl. (J2000) 

F225W 

F275W 

F336W 

F390W 

F435W 

F475W 

F555W 

F606W 

F625W 

F775W 

F814W 

F850LP 

F105W 

FllOW^ 

F125W 

F140W 

F160W 

IRAC chl 
IRAC ch2 



06:47:55.731 
+70:14:35.76 
-129 ± 51 nJy (-2.5cr) 
-95 ± 51 nJy (-1.9cr) 

2 ± 37 nJy ( O.Ocr) 
-8 ± 20 nJy (-0.4cr) 

± 26 nJy ( O.Ocr) 
-2 ± 14 nJy (-O.lcr) 
-3 ± 9 nJy (-0.3cr) 

3 ± 16 nJy ( 0.2cr) 
-35 ± 21 nJy (-1.7cr) 

4 ± 30 nJy ( 0.2cr) 
± 8 nJy ( O.lcr) 

-3 ± 30 nJy (-O.lcr) 
11 ± 12 nJy ( 0.9cr) 
-8 ± 10 nJy (-0.8cr) 
-3 ± 10 nJy (-0.3cr) 
63 ± 10 nJy ( 6.0cr) 
= 26.90 ±0.17 mag AB 
162 ± 13 nJy (12.4cr) 
= 25.88 ± 0.09 mag AB 

< 277 nJy^ 

< 245 nJy^ 



06:47:53. 
+70:14:22. 
-40 ± 50 nJy 
-31 ± 42 nJy 

49 ± 29 nJy 
1 ± 19 nJy 

43 ± 24 nJy 
-27 ± 16 nJy 

12 ± 7 nJy 

13 ± 20 nJy 
-52 ± 24 nJy 
-16 ± 52 nJy 

-2 ± 5 nJy 
1 ± 29 nJy 

14 ± 13 nJy 
3 ± 9 nJy 
7 ± 16 nJy 

50 ± 8 nJy 
= 27.15 ±0.17 
136 ± 9 nJy 
= 26.07 ±0.07 

< 166 
436 ± 139 nJy 



112 
94 
-0.8cr) 
-0.8cr) 
1.7cr) 
O.lcr) 
1.8cr) 
-1.7cr) 
1.7cr) 
0.6cr) 
-2.2cr) 
-0.3cr) 
-0.3cr) 
O.Ocr) 
l.lcr) 
0.3cr) 
0.5cr) 
6.7cr) 
mag AB 
(15.1cr) 
mag AB 
nJy 
( 3.1cr) 



06:47:55.452 
+70:15:38.09 



12 ± 32 nJy 



49 
-25 
1 
5 
7 
6 

-1 

23 
4 

-2 
6 
3 
7 
2 

26 



± 24 nJy 
± 18 nJy 
10 nJy 

14 nJy 
8 nJy 

4 nJy 
6 nJy 

10 nJy 
10 nJy 

3 nJy 

15 nJy 

5 nJy 

4 nJy 

5 nJy 
4 nJy 

3 ±0.17 mag AB 
4 nJy (lO.lcr) 
27.34 ±0.10 mag AB 

< 166 nJy 

< 138 nJy 



0.4cr) 
2.0cr) 
-1.4cr) 
O.lcr) 
0.4cr) 
0.9cr) 
1.4cr) 
-O.lcr) 
2.3cr) 
0.3cr) 
-0.8cr) 
0.4cr) 
0.6cr) 
1.9cr) 
0.4cr) 
6.1cr) 



= 27.^ 
42 ± 



± 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 

= 26.04 ± 
341 ± 16 
= 25.07 ± 

< 

436 ± 314 



-157 
-77 
25 
-6 
48 
-22 
15 
15 
-64 
-8 
-3 
4 
28 
2 
6 

139 



78 
70 
50 
29 
38 
23 
12 
26 
33 
61 
10 
45 
18 
14 
20 
14 



(-2.0cr) 
(-1.1a) 
( 0.5cr) 
(-0.2cr) 
( 1.3cr) 
(-l.Ocr) 
( 1.3cr) 
( 0.6cr) 
(-1.9cr) 
(-O.lcr) 
(-0.3cr) 
( O.lcr) 
( 1.6cr) 
( O.lcr) 
( 0.3cr) 
( 9.9cr) 
0.11 mag AB 
nJy (21.3cr) 
0.05 mag AB 
363 nJy 
nJy ( 1.4cr) 



nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 
nJy 



Note. — Fluxes in nanoJanskys (nJy) may be converted to AB magnitudes via uiab ~ 26 — 2.5 log2o(^i//(145 nJy)). Magnitude 
uncertainties, where given, are non-Gaussian but are approximated as 2.51og]^Q(e) times the fractional flux uncertainties. 

^ Sum of all three images with uncertainties added in quadrature. 

Visit A2 only, excluding visit A9 which exhibits signiflcantly elevated and non-Poissonian backgrounds due to Earthshine (§3.1). 
^ Includes uncertainties from modeling and subtracting a nearby brighter galaxy. More conservative estimates of these uncertainties 
were also considered in the analysis (§3.2). 



Table 3 

17 Clusters Searched in This Work 



High Magniflcation?^ 


Cluster^ 


Redshift 




AbeU 383 (0248.1-0331) 


0.187 




AbeU 611 (0800.9±3603) 


0.288 




AbeU 2261 (1722.5±3207) 


0.244 




MACSJ0329. 7-0211 


0.450 


Y 


MACSJ0647.8±7015 


0.591 


Y 


MACSJ0717.5±3745 


0.548 




MACSJ0744.9±3927 


0.686 




MACSJ1115.9±0129 


0.355 


Y 


MACSJ1149.6±2223 


0.544 




MACSJ1206. 2-0847 


0.439 




MACSJ1720.3±3536 


0.387 




MACSJ1931.8-2635 


0.352 


Y 


MACSJ2129.4-0741 


0.570 




MS2137.3-2353 


0.313 




RXJ1347.5-1145 


0.451 




RXJ1532.9±3021 


0.363 




RXJ2129.7±0005 


0.234 



^ CLASH clusters were selected based on either X-ray or strong 
lensing properties. The latter "high magniflcation" clusters are 
marked with Y's here. For details, see Postman et al. (2012). 
^ R.A. &L Decl. (J2000) are given in parentheses for the Abell clus- 
ters, encoded as they are in the names of the other clusters. 



filters. JDl is not detected above 1-a in any filter blue- 
ward of F140W. 

All three J-dropouts are confidently detected in two 
filters (F140W and F160W) observed at six different 
epochs over a period of 56 days (Fig. 5). No significant 
temporal variations are observed in position or bright- 
ness, ruling out solar system objects and transient phe- 
nomena such as supernovae, respectively (see Figs. 3 and 
6 and §5). 



3.1.3. Exclusion of F HOW Second Epoch 

Based on an initial standard reduction of the HST im- 
ages and standard SExtractor photometry, MACS0647- 
JD2 was detected in FllOW at 5-cr, while JDl and JD3 
were not significantly detected (0.9-cr and 1.7-cr, respec- 
tively). Our empirical rederivations of the photomet- 
ric uncertainties, including proper accounting for corre- 
lated noise (§3.1.1), reduced the significance of this detec- 
tion to 2.5-<j. However we ultimately we concluded this 
marginal detection was completely spurious due to sig- 
nificantly elevated and non-Poissonian backgrounds due 
to Earthshine in two out of five FllOW exposures, both 
obtained during the second epoch (see below). After ex- 
cluding these exposures, the detection significance drops 
to 0.3-cr, consistent with background noise. For refer- 
ence, see the WFC3/IR images in Fig. 2. 

Even based on the initial "standard" analysis described 
above, we determined that MACS0647-JD is at z < 9 
with a likelihood of '^10~^ based on a joint photomet- 
ric redshift analysis of all three images (§4). This like- 
lihood decreased further to 3xl0~^^ based on our im- 
proved analysis. These values are summarized in Table 
4. The spuriously high flux measurements may be seen 
in Fig. 3. 

The final observations of MACSJ0647.7+7015 were 
two 502-second exposures in FllOW obtained during 
visit A9, the second epoch for that filter. We found these 
to have significantly elevated backgrounds of 1.9 (sigma- 
clipped mean) ± 0.44 (RMS) and 6.4 ±0.27 electrons per 
second, respectively, compared to the more typical val- 
ues around 1.5 ±0.08. These high backgrounds were due 
to Earthshine, or sunlight reflected from the Earth. The 
first observation was obtained during twilight as the tele- 



6 



Coe et al. 2012 



Table 4 

Effects of FllOW Aberrant Second Epoch 



FllOW detection cr 


JDl 


JD2 


JD3 


P(z < 9)^ 


SExtractor photometry 


0.9 


5.0 


1.7 


1x10-9 


Empirical uncertainties 


0.5 


2.5 


1.0 


4x10-8 


Excluding second epoch 


-0.8 


0.3 


1.9 


3x10-13 



Note. — JD2 is spuriously detected in FllOW images processed 
using standard techniques. This is due to significantly elevated 
non-Poissonian backgrounds in the second epoch of observations 
due to Earthshine. We exclude this epoch in our analysis. See §3.1 

^ Based on the summed photometry of all three images, and as- 
suming MACS0647-JD is a galaxy well described by our templates. 
See §§4-5. 

scope pointed within 67-59 degrees of the bright hmb of 
the Earth. This Earth limb angle continued to steadily 
decrease from 47 to 24 degrees during the second obser- 
vation which was obtained during daylight. In the ob- 
servation log, the diagnostic Earth bright limb flag was 
raised halfway through the second exposure. We also ex- 
amined the ten individual readouts of 100 seconds each 
obtained over the course of both exposures and found 
the mean background increased steadily from 0.9 to 7.5 
electrons per second. The resulting elevated background 
RMS values of 0.44 and 0.27 electrons per second in the 
two exposures are the highest and sixth highest relative 
to the median values for a given filter in 1,582 CLASH 
observations to date of 17 clusters. None of the three 
F160W observations obtained at the beginning of visit 
A9 exhibit elevated backgrounds because they were ob- 
tained at night (twilight had yet to set in) and the Earth 
is less bright in F160W. 

Specifically, when we compared the measured RMS val- 
ues to what would be expected from scaling the back- 
ground intensity levels, we found that these RMS val- 
ues are several times higher than would be expected 
in the case of Poissonian statistics. We attribute this 
to the fact that the sky background was increasing 
in a strongly non-linear fashion during the exposure, 
whereas the up-the-ramp slope fitting algorithm imple- 
mented in "calwfS" implicitly assumes that the count 
rate is constant when converting measured counts into 
counts/second (see Dressel 2011). Since this assumption 
is violated, the the pixel-to-pixel variations in the final 
count-rate image no longer scale as expected for Poisso- 
nian statistics, as demonstrated by the much higher RMS 
values. Since these data no longer conform to Poissonian 
statistics, we were able to demonstrate that attempt- 
ing to combine them with the other data did not yield 
an improvement in S/N but instead produced combined 
datasets with non-Poissonian statistics, from which we 
were not able to obtain reliable photometry. 

We therefore exclude the two FllOW visit A9 expo- 
sures from our analysis and derive photometry instead 
from the weighted sum of the three visit A2 exposures. 

3.2. Spitzer Photometry 

To derive photometry in the longer wavelength Spitzer 
IRAC images (Fig. 7), we performed both GALFIT PSF 
fitting and aperture photometry on JD2 and JDS. No sig- 
nificant fiux is detected for either object in either chan- 
nel except for a 3-a detection of JD2 in ch2: mag = 
24.8 ± 0.3. Aperture photometry (2.4" diameter aper- 



ture) yields mag = 25.8 ±0.3, subject to an approximate 
0.7 mag correction, roughly consistent with the GALFIT- 
derived photometry. 

JDl is significantly contaminated by light from a 
nearby cluster galaxy. We modeled this galaxy using 
GALFIT, subtracted it from the image, and measured 
photometry in 2.4" diameter apertures, yielding a null 
detection plus uncertainty. We also added a simulated 
25th magnitude source and used GALFIT to derive its 
photometry. We conservatively combined the uncertain- 
ties from these two measurements in quadrature to yield 
total uncertainties (1-cr upper limits) of 277 and 245 nJy 
in chl and ch2, respectively {3-a limits of mag 24.2 and 
24.1). We also experimented with inflating these uncer- 
tainties further by one magnitude (3-<j limits of mag 23.2 
and 23.1). This would increase the JDl z < 9 likelihood 
(see §4) from -3x10"'^ to 2xl0"^ and the likelihood 
based on the integrated photometry of all three images 
from 3x10-^^ to 2x10"^ 

4. PHOTOMETRIC REDSHIFT 

We perform two independent analyses of the 
HST+Spitzer photometry to estimate the photometric 
redshift of MACS0647-JD. These two methods, BPZ 
(§4.1) and LePHARE (§4.2) were the top two perform- 
ers out of 17 methods tested in Hildebrandt et al. (2010). 
They yielded the most accurate redshifts with the fewest 
outliers given a photometric catalog for galaxies with 
known spectroscopic redshifts. 

According to our gravitational lensing models (§6), 
MACS0647-JD1, 2, and 3 are likely three multiple images 
of the same strongly lensed background galaxy. Thus, in 
this section, we present photometric redshift likelihoods 
for each individual image, as well as jointly for the two 
brighter images and for all three images. 

4.1. Bayesian Photometric Redshifts (BPZ) 

We used BPZ (Bayesian Photometric Redshifts; 
Benitez 2000; Coe et al. 2006) for our primary pho- 
tometric redshift analyses. We modeled the observed 
HST+Spitzer photometry of MACS0647- JD using model 
SEDs from PEGASE (Fioc & Rocca-Volmerange 1997) 
which have been significantly adjusted and recali- 
brated to match the observed photometry of galaxies 
with known spectroscopic redshifts from FIREWORKS 
(Wuyts et al. 2008). The FIREWORKS data set in- 
cludes 0.38-24/im photometry of galaxies down to mag 
~ 24.3 (5-cr i^-band) and spectroscopic redshifts out to 
z ^ 3.7. In analyses of large datasets with high quality 
spectra, this template set yields < 1% outliers, demon- 
strating that it encompasses the range of metallicities, 
extinctions, and star formation histories observed for 
the vast majority of real galaxies. (In §5.1 we explore 
a still broader range of galaxy properties using a syn- 
thetic template set which has not been recalibrated to 
match observed galaxy colors.) These templates include 
nebular emission lines as implemented by Fioc & Rocca- 
Volmerange (1997). 

The Bayesian analysis tempers the SED model quality 
of fit with an empirically derived prior P{z,T\m) on the 
galaxy redshift and type given its (delensed) magnitude. 
Our prior was constructed as in Benitez (2000) and up- 
dated based on likelihoods P{z,T\m) observed in COS- 
MOS (Ilbert et al. 2009), GOODS-MUSIC (Grazian et al. 
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Figure 2. The three images of MACS0647-JD as observed in various filters with HST. The leftmost panels show the summed 11-hour 
(17-orbit) exposures obtained in 8 filters spanning 0.4-0.9/im with the Advanced Camera for Surveys. The five middle columns show 
observations with the Wide Field Camera 3 IR channel in F105W, FllOW, F125W, F140W, and F160W, all shown with the same linear 
scale in electrons per second. The F125W images were obtained at a single roll angle, and a small region near JD2 was affected by- 
persistence due to a moderately bright star in our parallel observations immediately prior (see also Fig. 5). The right panels zoom in by a 
factor of 2 to show F110W+F140W+F160W color images scaled linearly between and 0.1 /j,Jy. 
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Figure 3. Flux measurements in the individual epochs observed over a period of 56 days. Filters are colored F160W (red), F140W 
(yellow), F125W (green), and FllOW (blue) as both individual data points and solid bands as determined for the summed observations. 
The FllOW exposures obtained in the second epoch (visit A9) were found to have significantly elevated and non-Poissonian backgrounds 
due to Earthshine (§3.1). These were excluded in our analysis; we adopted the FllOW fiuxes measured in the first epoch (visit A2). 



2006; Santini et al. 2009), and the UDF (Coe et al. 2006). 
According to this prior (extrapolated to higher redshifts) , 
all galaxy types of intrinsic (delensed) magnitude ~28.2 
are over 80 times less likely to be at z ~ 11 than z ^ 2. 
Thus our analysis is more conservative regarding high 
redshift candidates than an analysis which neglects to 
implement such a prior (implicitly assuming a flat prior 
in redshift). The prior likelihoods for MACS0647-JD are 
uncertain both due to the prior's extrapolation to z ^ 11 
and uncertainty in MACS0647-JD's intrinsic (delensed) 
magnitude. Yet it serves as a useful approximation which 
is surely more accurate than a flat prior. 



Based on this analysis, we derived photometric redshift 
likelihood distributions as plotted in Fig. 8 and summa- 
rized in Table 5. The images JDl, JD2, and JDS are 
best fit by a starburst SED at z ~ 10.9, 11.0, and 10.1, 
respectively. After applying the Bayesian prior, we find 
JDl and JD2 are most likely starbursts at 2: ~ 10.6 and 
11.0, respectively. A z ~ 2.5 elliptical template is slightly 
preferred for JDS, however z = 11 is within the 99% con- 
fidence limits (CL). Observed at mag ~ 27. S, we may not 
expect this fainter image to yield as reliable a photomet- 
ric redshift. 

In Table 4 we also provide joint likelihoods based on 
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Figure 4. Observed HST photometry (filled circles and triangles) plotted against the expected fiuxes (open blue squares) from a young 
starburst galaxy spectrum (gray line) redshifted to z ~ 11. HST filter transmission curves are plotted in the upper panel, normalized to 
their maxima, and with black dots indicating the effective "pivot" wavelengths. Photometry of the J-dropouts observed through these 
filters (Table 2) is plotted as the larger circles and triangles for positive and negative observed fiuxes, respectively, with 1-cr error bars. For 
some points, horizontal "error bars" are plotted to reiterate the filter widths. The gray line is a model spectrum of a young starburst at 
z = 11.0, the best fit to the summed photometry. The integrals of this spectrum through our filters give the model predicted fiuxes plotted 
as blue squares. Other galaxy types at 2; ~ 11 yield similar predicted HST fiuxes, as the shape of the spectrum cannot be constrained by 
the HST photometry alone. Redshifted Lyman-o; at 0.1216/im(l + z) ~ 1.46/im is indicated by the vertical dashed line. 



the brighter two images and ah three images equahy 
weighted. To properly downweight the fainter image, we 
also analyzed the integrated photometry of all three im- 
ages (with uncertainties added in quadrature). Based on 
this analysis including our Bayesian prior, and assuming 
MACS0647-JD is a galaxy well described by our template 
set (see also §5.1), we found z = 10.7^0.4 (95% CL) with 
a -3x10-^^ likelihood that MACS0647-JD is at z < 9. 
This likelihood corresponds to a 7.2-(j confidence that 



MACS0647-JD is at z > 9. The joint likelihood analysis 
(weighting all images equally) yields a similar 95% CL 
[10.2-11.1] and a more conservative P{z < 9) ~ 2xl0~^, 
or z > 9 at 5.5-<j. 

The strong confidence in the high redshift solution re- 
quires the combined HST and Spitzer photometry. With- 
out the Spitzer photometry, the z > 9 likelihood would 
drop to 95% for the summed HST photometry. Simi- 
larly, we would find P{z > 9) ~ 91% for JDl individ- 
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Figure 5. MACS0647-JD as observed in each of the individual epochs of F160W and F140W obtained over a 56-day period. These 
observations were obtained at two different telescope roll angles which alternate between the stamps shown here. A small region of the 
WFC3/IR images in our first roll angle was affected by persistence due to a moderately bright star in our parallel observations immediately 
prior. These pixels happen to fall within 1" of JD2 at that roll angle (marked in gray here and flagged as unreliable). Excluding this roll 
angle for JD2 does not signiflcantly affect the derived photometry. 
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Figure 6. Relative centroid measurements for the detections in F160W (red) and F140W (yellow) in individual epochs (circles) and 
summed observations (squares). Centroids measured in the summed NIR images are also plotted as gray diamonds. The offsets are generally 
less than one of our drizzled pixels (0.065^^), roughly half the native WFC3/IR pixel size (~0.13^^). 
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Figure 7. Spitzer IRAC chl (3.6/im) and ch2 (4.5/im) images of MACS0647-JD compared to the HST WFC3/IR F160W (1.6/im) image. 
Two intensity scalings and zooms are shown. Left: Both 26'' x 26'' F160W cutouts are scaled hnearly in photon counts to the same range 
as Fig. 2. And for each Spitzer filter, the same count range is used in each row. The background photon counts are significantly higher 
near JDl and JD2 (top row) due to intracluster light and scattered starlight. MACS0647-JD is not detected brightly in the Spitzer images, 
supporting the high-redshift solution. The only possible detection we report is for JD2 at 3.1-cr in ch2 (Table 2). JDl is contaminated by 
light from other nearby galaxies which we modeled and subtracted to estimate JDl's photometry. Right: In each of these 5" x 5" closeups, 
the intensity is scaled independently to the observed range within the central 3" x 3" . 
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Figure 8. Top row: photometric redshift probability distributions based on our BPZ analysis (§4.1) of HST+Spitzer photometry for 
each image. Cumulative probabilities P(< z) are shaded gray; probabilities P(z) per unit 0.01 in redshift are drawn as black lines; and 
likelihoods for individual SED templates are drawn as colored lines. In this work, we use four elliptical templates (Ell), one E/SO, two 
spirals (Sbc and Scd), and four starbursts (SB). We then interpolate nine templates between each pair of adjacent templates. Bottom left: 
joint likelihoods for all three images. Bottom center: likelihoods based on the summed photometry of all three images. For all of these 
likelihoods, we assume the prior plotted at bottom right for galaxies of intrinsic (delensed) magnitude ~28.2. This prior was empirically 
derived from large surveys with photometric and spectroscopic redshifts and extrapolated to higher redshifts. 
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Table 5 

Individual and Joint Redshift Likelihoods 



Image 95% CL 99% CL F(z < 9) 



multiple images according to our lens models (§6). Nev- 
ertheless, we found the only stars or brown dwarfs con- 
sistent with the observed colors are rare, transient post- 
AGB flare-ups, though these would be far more lumi- 
nous if observed within the Galaxy (§5.3). Solar system 
objects would have likely exhibited parallax motion and 
are inconsistent with the observed colors (§5.4). Interme- 
diate redshift long duration multiply-imaged supernovae 
(§5.5) and emission line galaxies (§5.6) are also extremely 
unlikely. We conclude that MACS0647-JD is most likely 
either at z ~ 11 or exhibits unique photometry yet to be 
observed in any other known object. 

5.1. Intermediate Redshift Galaxy? 
5.1.1. SED Constraints 

While we found P{z < 9) - 3x10"^^, the next best 
alternative to z 11 is an early type galaxy (ETC) 
at z - 2.5 (Fig. 9). At ^ - 2.65, the 4000A break is 
redshifted to 1.46/im, coinciding with Lyman-a (12 16 A) 
redshifted to z ~ 11.0. However, 4000 A breaks are not 
expected to be as strong as observed for MACS0647-JD 
(Figs. 10 and 11). JDl features a J125 — i^ieo ^ 3 magni- 
tude break between F125W and F160W as well as a ~1 
magnitude break between F140W and F160W. Thus low 
redshift ETGs yield a significantly worse SED fit than 
z ~ 11 for all three images as quantified in Fig. 8 and 
Table 5. 

To explore an even broader range of galaxy SED mod- 
els than used in §4, we utilized the flexible stellar popula- 
tion synthesis (FSPS) models from Conroy et al. (2009) 
and Conroy & Gunn (2010). They provide simple stel- 
lar population (SSP) models which span ages of 5.5 < 
log(age/yr) < 10.175 and metallicities of 0.0002 < Z < 
0.03 (where Zq = 0.019). Nebular emission lines are 
not included. We convolved their SSP models with star 
formation histories (SFH) ranging from the single early 
burst (SSP) to exponentially declining ("r models"), con- 
tinuous (constant rate), and exponentially rising ("in- 
verted r models"). The latter rising SFH likely describes 
high-redshift galaxies best according to both observa- 
tions (Maraston et al. 2010; Papovich et al. 2011; Reddy 
et al. 2012) and simulations (Finlator et al. 2011). Fi- 
nally we added a variable degree (up to Ay = 30 mag- 
nitudes) of Calzetti et al. (2000) dust extinction with 
Rv = Av/Eb-v = 4.05. 

To uncover the most likely solutions in different re- 
gions of this multi-dimensional parameter space, we be- 
gan with relatively coarse grid searches with redshift in- 
tervals of 0.1 and ~9 steps in each of the four other free 
parameters. We then zoomed in on the higher likelihood 
regions, found again to be roughly z ^ 2.5 and 11. Fi- 
nally, we ran Powell (1964) minimizations to find the best 
fitting SEDs at each of these redshifts. 

We supplemented these SEDs with a suite of smooth 
r models with stochastic bursts superposed (e.g., Kauff- 
mann et al. 2003; Salim et al. 2007), as well as truncated 
("quenched") SFHs designed to reproduce the colors of 
post-starburst (K+A) galaxies. 

Our results with this combined template set confirm 
that a z - 11 model fits MACS0647-JD best, while 
evolved and/or dusty galaxies at z ~ 2.5 provide the 
best alternatives but are still significantly worse statis- 
tically. The best fitting intermediate-redshift template 



JDl (F160W ~ 25.9)^ 10 

JD2 (F160W ~ 26.1)^ 10 

JD3 (F160W ~ 27.3)^ 2 

P(JDl)xP(JD2)b 10 
P(JDl) xP(JD2) xP(JD3)b 10 

P(JD1 + JD2)^ 10 

P(JD1 + JD2 + JD3)^ 10 



.62+g-J^ [10.11-11.67] 3x10-'^ 

•99lg;^° [9.99-11.69] 3x10"^ 

AStoll [1.81-11.07] 7x10-1 

.66lg;^^ [10.21-11.53] 4x10-^0 

•42+g;^^ [10.12-11.30] 2x10"^ 

•99+g;^^ [10.22-11.59] 2x10'^^ 

.7llo;37 [10.20-11.47] 3x10-^3 



Note. — BPZ results assuming MACS0647-JD is well modeled 
by our SED templates (§§4-5). All likelihoods include a Bayesian 
prior which assumes galaxies of this (unlensed) magnitude are over 
80 times more likely to be at z ~ 2 than 2; ~ 11. See also Fig. 8. 

^ Approximate AB magnitudes are given in parentheses. (See also 

Table 2.) Note JD3 is significantly fainter. 

^ Joint likelihood of multiple images weighted equally. 

^ Likelihood based on integrated photometry of multiple images. 

ually. However, the most likely solutions for JD2 and 
JDS would be early types at z ~ 4. We would expect 
such galaxies to be mag ~ 23 in the Spitzer observa- 
tions, which is extremely unlikely (as quantified above) 
given the measured photometry (see also §3.2). 

4.2. Le PHARE 

We also used Le PHARE (Arnouts et al. 1999; Ilbert 
et al. 2006, 2009) to independently estimate the photo- 
metric redshifts. For this analysis we used a SED tem- 
plate library primarily from Ilbert et al. (2009) as opti- 
mized for the COSMOS survey (Scoville et al. 2007b,a; 
Koekemoer et al. 2007). This includes three ellipticals 
and seven spirals as generated by Polletta et al. (2007) 
using the GRASIL code (Silva et al. 1998), as well as 12 
starburst galaxies with ages ranging from 30 Myr to 3 
Gyr generated by GALAXEV based on Bruzual & Char- 
lot (2003). We supplemented these with four additional 
elliptical templates for a total of seven ellipticals. 

We added dust extinction in ten steps up to E{B — 
V) = 0.6. (Stronger degrees of extinction are explored in 
§5.1.1.) Four different dust laws were explored: Calzetti 
et al. (2000); Calzetti plus two variations on a 2 170 A 
bump; and Prevot et al. (1984) as observed for the SMC. 

We adopted the Benitez (2000) prior as implemented 
in LePhare. The results were consistent with those from 
BPZ: z = lOStoi (JDl), z = 10.6t^-j (JD2), and z = 

10- 1^0.3 (JD3), each at 68% CL. A secondary solution of 
z ^ 2.5 was reported for JD3 with a peak likelihood ten 
times less than that of the best fit high redshift solution. 

5. LOWER REDSHIFT INTERLOPERS RULED OUT 

In this section we consider a broad range of 2; < 11 
possibilities. As found in §4, the z < 9 likelihood is 
formally ~3xl0~^^ assuming MACS0647-JD is a galaxy 
well-modeled by our SED templates. Though strongly 
disfavored, a z ~ 2.5 early type and/or dusty galaxy 
is the most likely alternative, as we discuss further in 
§5.1. We reanalyzed previously published J-dropouts 
and found them most likely to be at intermediate red- 
shift (§5.2). Objects within the Galaxy are less likely, as 
this would require three objects with extremely rare col- 
ors (Fig. 13) at positions consistent with strongly- lensed 
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Wavelength A (^m) Wavelength A (//m) Wavelength A (^m) 

Figure 9. Observed near-IR photometry from HST WFC3/IR and Spitzer IRAC (filled circles and triangles) compared to the expected 
fiuxes (open squares) from two SEDs: the z = 11.0 starburst from Fig. 4 (blue) and a z = 2.5 early type galaxy (red). Note the JD3 plot 
is scaled differently along the y-axis. 
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Figure 10. Observed WFC3/IR colors (shaded 68% confidence regions) for JDl (solid lines), JD2 (dashed lines), and JDS (dotted lines) 
plotted against those predicted with the BPZ template set from young starburst (blue) to early type (yellow-orange-red) as a function 
of redshift. The three panels plot fiux ratios in F125W / F160W (left), F125W / F140W (middle), and F140W / F160W (right). The 
corresponding colors in magnitudes are given along the right axes. 

The uncertainties on the z ^ 11 SED parameters are 
quantified in §7.3. A proper calculation of the red- 
shift likelihoods based on these templates would require 
an estimate of the prior likelihoods in this multidimen- 
sional parameter space, which is beyond the scope of this 
work. And while these templates probe a broad param- 
eter space, we derive our primary photometric redshift 
estimates in §4 from templates which have been well cal- 
ibrated to match the observed photometry of galaxies 
with spectroscopic redshifts. 

We note there is no evidence that z > 2 ETGs have 
significantly different SEDs than our ETG models cali- 
brated at lower redshifts. The highest redshift ETG ob- 
served to date is HUDF-1446 with a spectroscopic red- 
shift z = 2.67 (Damjanov et al. 2011). Coe et al. (2006) 
published a photometric redshift z = 2.74 ± 0.44 for this 
object using BPZ, in good agreement with the true red- 
shift. Their ETG templates yielded a good fit to the 
ACS (^435^606^775450) ^ud NICMOS {JnoHieo) pho- 
tometry, including the Jno — i^i60 = 1.89 ± 0.13 break 
with HiQo = 23.074 ± 0.098 and significant detections in 
ah ACS filters. 




-2 -1 01234567 
Color (magnitudes): F160W - fRACl 

Figure 11. Observed colors in WFC3/IR F125W - F160W and 
F160W - IRAC chl plotted as black lines (95% confidence limits) 
versus those predicted from the current BPZ template library (lines 
colored as a function of redshift and made thicker for earlier galaxy- 
types). 



to the summed photometry {z ~ 2.7; ~400 Myr old; 
Ay ~ 0.8 mag) with = 57.6 is only ~10~^ times as 
likely as the best fitting z ^ 11 template {z ~ 10.9; ^6 
Myr old; Ay = 0) yielding = 16.9 with > 14 de- 
grees of freedom given the 19 photometric measurements 
and < 5 free parameters (see discussion in Andrae et al. 
2010). 



5.1.2. Lower Stellar Mass Than Observed z >2 ETGs 

If MACS0647-JD were at z - 2.5 (despite the low like- 
lihood of this from SED fitting) it would likely be the 
least massive early type host galaxy observed to date at 
z > 1. Spectroscopically confirmed z > lA ETGs to date 
have stellar masses > 2xlO^^M0 (Damjanov et al. 2011). 
HUDF-1446 at z = 2.67, for example, is -SxlO^^M©. 
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Our subset of lens models that allow for MACS0647- 
JD to be at z ~ 3 (§6.3) suggest that the magnification 
of the brightest two images would be /i > 30. Thus it 
would be intrinsically ~300 times fainter than HUDF- 
1446, with a correspondingly lower stellar mass on the 
order of ~2xlO^M0 (and stih < lO^M© if we assume 
a more conservative magnification factor of /i ~ 10; see 
also z ^ 11 mass estimates in §7.2). 

Quiescent galaxies of such low masses at z > 2 would 
be a surprising discovery. Observations to date demon- 
strate (e.g., Peng et al. 2010) that star formation is 
only significantly quenched by feedback in more massive 
galaxies, or alternatively as a galaxy is harassed as a 
satellite of a larger halo. MACS0647-JD is not observed 
to be a satellite of a galaxy group. 

5.2. Comparisons to Previously Published J -dropouts 

The previous highest redshift candidate, UDFj- 
39546284 (Bouwens et al. 2011a), was detected at 5.8-cr 
in a single HST band (WFC3/IR F160W) dropping out 
of F125W and bluer filters also with non-detections in 
Spitzer yielding a photometric redshift oi z = 10.3 ±0.8. 
The ultimate inclusion of the F140W filter on WFC3 
(Brown & Baggett 2006) and in the CLASH observing 
program enable us to securely identify MACS0647-JD 
as the highest redshift galaxy candidate to date. At 
z ^ 11.0, Lyman-a is redshifted to ~1.46/im, causing 
the galaxy light to drop out of ^2/3 of the F140W band- 
pass as well as ~l/5 of F160W. The ratio between these 
two filling factors (0.8/0.33 ~ 2.4, corresponding to ~1.0 
mag) places tight, model-independent constraints on the 
wavelength of the (redshifted) Lyman break and thus the 
redshift of MACS0647-JD (Fig. 10). The five NIR HST 
filters used by CLASH also enabled Zheng et al. (2012) 
to discover a J-dropout lensed by M ACS Jl 149.6+2233 
and robustly measure its photometric redshift to be 
z = 9.6 ±0.2 (68% CL). 

Laporte et al. (2011) identified a J-dropout lensed 
by Abeh 2667 based on VLT (F0RS2 and HAWK-I), 
ACS/F850LP, and Spitzer IRAC (chl through ch4) pho- 
tometry. Hayes et al. (2012) then measured a spectro- 
scopic redshift of z = 2.082 for that galaxy, A2667-J1. 
Laporte et al. (2011) had already emphasized that z > 9 
possibilities were excluded based on the significant (6.0- 
cr) ACS detection. We concur with this conclusion af- 
ter reanalyzing their photometry as provided in Hayes 
et al. (2012). Only by excluding the ACS data point 
and assuming no Bayesian redshift prior do z > 9 so- 
lutions have significant probability (Fig. 12). If Spitzer 
IRAC ch3 and ch4 were not available (as is the case with 
MACS0647-JD) in addition to the ACS detection being 
unavailable, then the z > 9 likelihood would rise further 
yet still be insignificant once the prior is included. The 
z > 9 likelihood is enhanced further, but only modestly, if 
the IRAC chl and ch2 uncertainties are infiated to yield 
only 3-<j detections (as is the case for our JD2 IRAC 
ch2). In Fig. 12 we compare the observed NIR photom- 
etry of A2667-J1 and MACS0647-JD1. Our multiband 
HST photometry of the latter yields significantly tighter 
upper limits on the non-detections and adds a key data 
point at 1.4/im, resulting in a far greater z ^ 11 likeli- 
hood even when accounting for the Bayesian prior which 
disfavors them (Fig. 8). 

We also applied our analysis methods to the photom- 



etry of other J-dropouts in the literature. Schaerer 
et al. (2007) showed A1835-#17 was fit weh by a dusty 
{Ay ~ 3.6 mag) starburst at z ~ 0.8. Dickinson et al. 
(2000) presented both z > 2 and z > 10 solutions for 
HDF-N J123656.3±621322. And HUDF-JD2 (Mobasher 
et al. 2005) has since been shown to likely be a 2; ^ 1.7 
LIRG (Chary et al. 2007). For ah three of these J- 
dropouts, our analysis yields low redshift (2 < z < 4 or 
very dusty ^ ^ 1) solutions which are strongly preferred 
given our Bayesian prior. 

5.3. Stars or Brown Dwarfs? 

MACS0647-JD1, JD2, and JD3 are most likely multi- 
ple images of a strongly lensed background galaxy, well 
behind the z = 0.591 cluster. Their observed colors are 
extremely rare in our multiband HST catalogs of 17 clus- 
ters observed to date (Fig. 13). And they lie at or near 
the predicted positions of multiply- lensed images (§6). 
It would be highly unlikely to find three foreground (un- 
lensed) objects with such rare colors coincidentally at 
these positions. Still we consider here possible interlop- 
ers within the Galaxy, namely stars, brown dwarfs, and 
(in §5.4) solar system objects including KBOs (Kuiper 
Belt Objects) and Oort Cloud objects. 

JDl and JD2 are perhaps marginally resolved with 
deconvolved FWHM < 0.2^' (0.3^^ observed with a 0.2'^ 
PSF). We performed two independent analyses attempt- 
ing to determine whether the observed FWHM was large 
enough to definitively distinguish it from the stellar lo- 
cus. These analyses reached different conclusions. There- 
fore we turn to other lines of evidence to rule out stars 
and smaller objects. 

Stars are relatively plentiful in this field as the Galactic 
latitude is relatively low (±25.1°). We used the online 
tool TRILEGAL^i (Girardi et al. 2005) to calculate that 
we may expect ~5 late type M dwarfs of ~26th magni- 
tude or fainter within our FOV. However the predicted 
colors are J125 — i^ieo ~ 0.4, a break significantly weaker 
than that observed. 

In Fig. 14, observed and expected colors of stars and 
brown dwarfs (including types M, L, T, and Y) are plot- 
ted versus those observed for the J-dropouts. No dwarf 
color is able to reproduce the observed J-dropout col- 
ors. According to models, the colors of extremely cold 
(~200K) Y dwarfs come close to matching the red ob- 
served HST NIR colors, but these are expected to be 
significantly brighter in IRAC by up to 10 magnitudes. 
The coldest dwarfs yet discovered are Y dwarfs includ- 
ing WISEP J1828±2650 at -300 K (Gushing et al. 2011; 
Kirkpatrick et al. 2012) with colors as plotted in Fig. 14 
from ground-based JH and WISE W2 4.6/im observa- 
tions. 

Of the stellar spectra observed with IRTF (Gushing 
et al. 2005; Rayner et al. 2009), the M8III red giant WX 
Piscium (IRAS 01037±1219; Ulrich et al. 1966; Decin 
et al. 2007) comes closest to matching the observed col- 
ors of MACS0647-JD. However such a large, bright star 
(M — —4) would need to be well outside the Galaxy 
(— 10 Mpc distant) to be observed at 26th magnitude in 
F160W (as argued in Dickinson et al. 2000 and Bouwens 
et al. 2011a for two previous z ^ 10 candidates). If 
MACS0647-JD were within the Galaxy (out to -10 kpc), 

http://stev.oapd.inaf.it/trilegal 
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Figure 12. Left: comparison of near infrared photometry of lensed J-dropouts MACS0647-JD1 (this work) and A2667-J1 (Laporte 
et al. 2011; Hayes et al. 2012). Also overplotted are three SED fits from Figure 2 of Hayes et al. (2012) to the photometry of A2667-J1 
at its spectroscopic redshift z = 2.082. The Hayes et al. (2012) photometry plotted here is all from VLT (FORS2 and HAWK-I) except 
for the 6.0-cr detection in HST/ACS F850LP, the upper diamond with a darker border at 0.91^im. Right: photometric redshift probability 
distribution for A2667-J1 based on our reanalysis of the photometry provided in Hayes et al. (2012) with and without a Bayesian prior. 
The top panel uses all the available photometry. The middle panel omits the 6.0-cr detection in ACS/F850LP. The bottom panel omits 
both ACS and Spitzer IRAC ch3 and ch4. The spectroscopic redshift z = 2.082 is indicated by the red vertical lines. 
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Figure 13. Observed NIR colors of the J-dropouts (red dia- 
monds with 1-cr uncertainties) plotted against those observed for 
all other 20,746 CLASH sources brighter than 28th magnitude in 
both F160W and F140W and also observed in F125W (filled circles 
and density map). The horizontal axis gives the ratio of F160W 
fiux to the maximum fiux in all bluer WFC3/IR filters. The verti- 
cal axis gives a similar fiux ratio but for F140W. Three objects with 
colors similar to the J-dropouts appear to be spurious IR artifacts 
based on visual inspection, and we mark these with X's. 

it would have an absolute magnitude of M +11 or 
fainter, consistent with a red dwarf in terms of magni- 
tude but not color as shown above. 

A few red giants in the post-AGB phase have been ob- 
served to flare up apparently as the result of a helium 
burning "thermal pulse" which triggers the ejection of a 
dust shell. "Sakurai's object" (Duerbeck & Benetti 1996) 
and WISE J1810-3305 (Gandhi et al. 2012) do have sim- 
ilar colors to our J-dropouts. But again these are very 
bright events, observed at 0.34 ±0.01 Jy and 2.74 ±0.06 
Jy, respectively, in the JJ-band with 2MASS. Believed 
to be a few kpc distant, they would need to be removed 
to several Mpc to be observed at ^0.1 fiJy as our J- 
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Figure 14. Observed colors in J125, Hiqq, and at 4.5/im plotted 
as black lines (95% confidence contours) versus those observed and 
predicted for stars and brown dwarfs. Colors derived from stellar 
spectra observed with IRTF (Gushing et al. 2005; Rayner et al. 
2009) are plotted as open magenta circles for dwarfs and open 
black star symbols for giants and supergiants. Blue error bars 
are observed photometry (ground-based and WISE) for Y dwarfs 
(Kirkpatrick et al. 2012). Open black diamonds are post-AGB fiare 
ups with dust ejecta observed with 2MASS and WISE; the upper 
diamond is "Sakurai's object" (Duerbeck &; Benetti 1996) and the 
lower diamond is WISE J1810-3305 (Gandhi et al. 2012). Simu- 
lated dwarf spectra from Hubeny & Burrows (2007) are plotted as 
filled circles colored as a function of temperature. 

dropouts. These are also rare events, lasting on the or- 
der of 100 years (but varying more rapidly), such that 
only these two have been reported to date. It would be 
highly unlikely to detect three such events occurring at 
the same time in the same HST field. 

5.4. Solar System Objects? 

If the J-dropouts were solar system objects, we would 
expect to have detected their proper motions in our 
6 epochs of F160W/F140W imaging spanning 56 days 
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(Figs. 5 and 6). Only an Oort cloud object at ~50,000 
AU would be orbiting the Sun sufficiently slowly for us 
not to have detected its motion. But Oort cloud objects 
are expected to be significantly fainter (~58th magnitude 
assuming a diameter of ~20 km; Sheppard 2010) and 
have different colors than those observed here. Benecchi 
et al. (2011) measured HST/NICMOS F110W-F160W 
colors of 80 Trans-Neptunian Objects and found they 
have HST F110W-F160W colors clustered around -0.6 
with none redder than 0.8. For Oort cloud objects to be 
observed at ~26th magnitude due to refiected sunlight 
alone, even with 100% albedo, their radii would have to 
be ~10^^ m which is larger than the Sun. Oort cloud 
brown dwarfs emitting thermally would be far brighter 
than observed and, as discussed in §5.3, have different 
colors. 

5.5. Lensed Supernova? 

The J-dropouts do not exhibit any significant tempo- 
ral variations in brightness either over our 56 days of 
observations (Fig. 3), ruling out most transient phenom- 
ena such as supernovae. However, Type IIP supernovae 
can plateau to a roughly constant magnitude for ~100 
days (e.g., Arcavi et al. 2012) which we would observe 
to last ~100(1 + z) days due to cosmic time dilation 
(Blondin et al. 2008). A Type la supernova at z ~ 4 
would be observed to have magnitudes and colors simi- 
lar to those observed in HST for MACS0647-JD. A Type 
IIP plateau supernova would likely be bluer, but could 
perhaps match the observed HST colors. We would ex- 
pect to detect it as a bright object at longer wavelengths, 
but the Spitzer images were obtained 1.5 years earlier, 
perhaps before the star went supernova. 

This intriguing scenario is ruled out by the gravi- 
tational lens time delays due to MACSJ0647.7+7015, 
which we estimate to be on the order of 1-10 years be- 
tween JDl and JD2 and ~50 years between these and 
JD3. Our subset of lens models (§6) which allow for 
z ^ 4 also suggest that the intrinsic fiuxes of all three 
images are roughly consistent with one another (at least 
to within a magnitude). Thus the supernova plateau 
(several magnitudes brighter than the host galaxy) would 
have to have lasted 50/(1+2:) ~ 10 years for us to observe 
it simultaneously in all three images. Even if this were 
somehow possible, we would then expect to have detected 
the earlier (least time-delayed) images with Spitzer. 

5.6. Emission Line Galaxy? 

In principle, an AGN / starburst galaxy with an un- 
detected continuum and two or more extremely strong 
nebular emission lines redshifted into F140W and F160W 
could reproduce the observed HST colors. The only plau- 
sible configuration is that H/3 (4861 A) and [OIII] (4959A) 
are redshifted to within F140W and F160W, while [OIII] 
(5007A) is redshifted beyond F140W but within F160W. 
This is possible for the narrow redshift range 2.20 < z < 
2.22. At this redshift, F140W and F160W have rest- 
frame widths of ~ 1229 A and '^889 A, respectively. In 
the case of JDl, we measure the fiux blueward of F140W 
to be —3.4 ±3.7 nJy. We conservatively adopt < 7.4 nJy 
as the 2-cr upper limit on the continuum fiux. Boost- 
ing the F140W flux to the observed ~63 nJy would re- 
quire emission lines with a combined equivalent width 



EW > 1229 X (63/7.4-1) - 9525 A. Similarly, increasing 
the F160W flux to the observed '^162 nJy, would require 
a combined EW > 889 x (162/7.4- 1) - 19114A (-10"^^ 
erg/s/cm^). Thus in our conflguration assuming a con- 
tinuum flux of 7.4 nJy: 

EW(H/3 + [OIII]4959) - 9234A, (1) 

EW(H/3 + [OIII]4959 + [OIII]5007) ^ 18573A, (2) 

EW( [OIII] 5007) -{2)-{l)^ 9339A, (3) 

EW([OIII]4959) ^ (3) / 3 ^ 3113A, (4) 

EW(H/3) = (1) - (4) ^ 6121A, (5) 

where the line ratio in Equation 4 is dictated by the 
relative transition probabilities. 

An [OIII] (5007A) line with EW > 9000A would be 
several times greater than the strongest emission lines 
observed to date, approaching EW ~ 2000A, for [OIII] 
(5007A) and Ra (6563A) (Atek et al. 2011; van der Wei 
et al. 2011; Shim et al. 2011; Fumagalh et al. 2012). If 
we consider instead the 5-cr continuum limit of 18.4 nJy, 
we would stih require EW - 3959A for [OIII] (5007A) 
and EW ~ 1659A for H/3. The strongest H/3 lines are ro- 
bustly predicted to have EW < 800A even for extremely 
young stellar populations, according to models like Star- 
burst9922 (Leitherer et al. 1999, 2010). Finally, given 
such bright lines, one would also expect a signiflcant con- 
tribution of [OH] (3727A) in FllOW and F125W, which 
is not observed. 

6. GRAVITATIONAL LENS MODELING 

We identifled 24 strongly-lensed images of 9 back- 
ground galaxies (§6.1), used them to derive lens models 
using three different methods (§6.2), and derived results 
including magniflcations in (§6.3). Importantly, our lens 
models show that MACS0647-JD1, 2, and 3 are observed 
in relative positions as expected if they are strongly- 
lensed multiple images of the same galaxy at z ~ 11. 

6.1. Strongly Lensed Multiple Images 

Zitrin et al. (2011a) presented a preliminary gravita- 
tional lens mass model of the MACSJ0647.7+7015 
cluster core based on pre-CLASH HST/ ACS 
F555W+F814W imaging and their identiflcations 
of two background galaxies strongly lensed to produce 
multiple images. Based on CLASH imaging in 15 ad- 
ditional HST fllters and additional lens modeling using 
the Zitrin et al. (2009) method, we have now identifled 
7 more galaxies which have been multiply imaged, and 
we have measured robust photometric redshifts for all 9 
galaxies. This enables us to model the mass distribution 
(primarily dark matter) and thus lensing properties in 
greater detail. 

In addition to the three images of MACS0647-JD at 
z ^ 11, we observe 21 multiple images of 8 back- 
ground galaxies with photometric redshifts ranging from 
2 < ^ < 6.5 (Table 6 and Fig. 1). The candidate z - 6.5 
system is notable in its own right, consisting of two im- 
ages observed at magnitudes ~ 26.3 and 27.3 in the NIR. 

http: / /www. stsci.edu /science /starburst 99 / 
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For each of systems 3, 5, and 8, our lens models pre- 
dict a third faint counterimage, but we cannot unam- 
biguously identify it among several possible candidates. 
To be conservative, we do not include these uncertain 
identifications as constraints on our lens models. Inclu- 
sion of these candidates does not significantly affect the 
lens models. In Fig. 1 we indicate the predicted locations 
of these counter images. 

6.2. Strong Lens Modeling Methods 

Based on the observed positions of all 24 
strongly lensed images, we model the mass of 
MACSJ0647.7+7015 using three different methods: 
the Zitrin et al. (2009) method, Lenstool (Kneib et al. 
1993; Julio et al. 2007), and LensPerfect (Coe et al. 
2008, 2010). The first two methods are "parametric" 
in that they assume light traces mass, which has 
proven to be a very good prior. For example, some 
of the earliest efforts to model cluster lenses found 
that assigning masses to individual luminous cluster 
galaxies significantly improved the reproduction of 
strongly lensed images (Kassiola et al. 1992; Kneib et al. 
1996). LensPerfect makes no assumptions about light 
tracing mass, exploring a broader range of mass models 
and perfectly reproducing the observed positions of all 
multiple images positions as input. 

The Zitrin et al. (2009) mass model parameteriza- 
tion consists of three components: the cluster galaxies, 
a smooth cluster halo, and an external shear. Cluster 
galaxies were identified according to the "red sequence" 
in F814W-F555W color-magnitude space, then verified 
with photometric redshifts. Each cluster galaxy was 
modeled as a power law density profile, its mass scal- 
ing with fiux observed in F814W. In this work we also 
allowed the masses of the two brightest central galax- 
ies to vary independently. The cluster halo component 
was derived from this galaxy component by smoothing 
the latter with either a polynomial spline or a Gaussian. 
The two components were allowed to scale independently 
before being added. In all, there were eight free param- 
eters: the mass scalings of the galaxy and halo compo- 
nents, the masses of the two brightest central galaxies, 
the power law of the galaxy density profiles, the degree 
of the polynomial spline or Gaussian smoothing width, 
and the amplitude and direction of the external shear. 

The Lenstool model consisted of an ellipsoidal NFW 
halo (Navarro et al. 1996) plus cluster galaxies modeled 
as truncated PIEMDs (pseudo-isothermal elliptical mass 
distributions; Kassiola & Kovner 1993). We assumed 
core radii Tcore = 300 pc and luminosity scaling relations 
as in Julio et al. (2007): velocity dispersion ctq cx L^/^ and 
cutoff radius rcut L^/^, resulting in all galaxies having 
equal mass-to-light ratios. The normalizations of these 
two scaling relations were free parameters along with the 
cluster halo position (x, ?/), ellipticity (e, ^), scale radius, 
and concentration. There were eight free parameters in 
ah. 

6.3. Lens Model Results 

Given the observed position of any one of the 
MACS0647-JD images, the Lenstool model accurately 
predicts and reproduces the positions of the other two im- 
ages to an RMS of 1.3'^ as minimized for z = 11.59l5;53. 



Table 6 

24 Multiple Images of 9 Strongly Lensed Galaxies 







R.A. 


Decl. 


Photometric 




ID 


( 72000 0"! 


( T2000 0"! 






la 


06 47 51.87 


+70 15 20.9 


2.2 ± 0.1 




b 


06 47 48.54 


+70 14 23.9 


2.2 + 0.1 




c 


n^j AT ceo ni 
Ub 47 oz.Ul 


1 Tr\ 1/1 C^Q o 

+70 14 oo.o 


2.2 i 0.1 




2a 


06 48 00.33 


+70 15 00.7 


4.7 + 0.1 




b 


06 48 00.33 


+70 14 55.4 


4.7 + 0.1 




c 


r\Ci AT C^Q /^O 

Ob 47 oo.bz 


1 Tr\ 1/1 oi o 
+70 14 zl.o 


4.7 + 0.1 




3a 


06 47 53.85 


+70 14 36.2 


3.1 + 0.1 




b 


06 47 53.41 


+70 14 33.5 


3.1 + 0.1 




4a 


06 47 42.75 


+70 14 57.7 


1.9 ± 0.1 




b 


06 47 42.93 


+70 14 44.5 


1.9 + 0.1 




c 


r\Ci AT /IK QT 

Ob 4 / 4o.o / 


1 Tr\ 1 K OK Q 

+ / io zO.o 


1 n 1 n 1 
i.y zL O.i 




5a 


06 47 41.04 


+70 15 05.5 


6.5 + 0.15 




b 


06 47 41.16 


+70 14 34.4 


6.5 + 0.15 


J JJi 


oa 


r\Ci AT ce ce TA 
Ob 4 00. 1 4 


1 Tr\ 1/1 Q K 7 
+ < 14 oO. 1 


11 n 1 n 9 
ii.O zL O.o 


JD2 


b 


06 47 53.11 


+70 14 22.8 


11.0 + 0.3 


JD3 


c 


06 47 55.45 


+70 15 38.0 


11.0 + 0.3 




7a 


06 47 50.91 


+70 15 19.9 


2.2 + 0.15 




b 


06 47 47.73 


+70 14 23.2 


2.2 + 0.15 




c 


06 47 48.69 


+70 14 59.8 


2.2 + 0.15 




8a 


06 47 48.61 


+70 15 15.8 


2.3 + 0.1 




b 


06 47 47.34 


+70 15 12.5 


2.3 + 0.1 




9a 


06 47 43.79 


+70 15 00.4 


5.9 + 0.15 




b 


06 47 44.98 


+70 14 23.2 


5.9 + 0.15 




c 


06 47 49.06 


+70 15 37.7 


5.9 + 0.15 



^ Joint likelihoods for each system with approximate 68% 
uncertainties. 



This scatter is consistent with the ^lA" expected due 
to lensing by hne of sight structures and variation in the 
mass-to-hght ratio of cluster galaxies (Julio et al. 2010; 
Host 2012). 

The lens model and inferred redshift for MACS0647-JD 
do not change significantly if the MACS0647-JD images 
are excluded as constraints. In this case, the 8-parameter 
lens model remains well constrained by the 21 other mul- 
tiple images which provide 26 constraints (see discussion 
below). 

Using the Zitrin et al. (2009) method, two sets of ac- 
ceptable models are found in different regions of the 
model parameter space. One set prefers 2; ^ 11 for 
MACS0647-JD, while the other prefers z - 3.5. The 
latter mass models have fiatter profiles. 

Our LensPerfect analysis confirms this is a degeneracy 
between the MACS0647-JD redshift and the cluster mass 
distribution. A wide range of redshifts including z = 3.5, 
11.0, and 11.6 is permitted given the strong lensing data. 
When fixing the redshift to any of these values, LensPer- 
fect produces reasonable lens models (physical and with 
light approximately tracing mass) which perfectly repro- 
duce all 24 observed positions of the 9 strongly lensed 
galaxies. When the MACS0647-JD redshift is set lower, 
the cluster mass distribution is more spread out yielding 
a fiatter profile. We confirm that the parametric models 
have similar differences, in part due to their parameteri- 
zations of the cluster mass distribution. 

Including the redshift of MACS0647-JD, both para- 
metric models have < 9 free parameters. (This num- 
ber should be considered a maximum given covariances 
among the parameters. See discussion in Andrae et al. 
2010.) There are 30 constraints = 2 x (24 - 9), where 
the constraints are the two coordinates (x^y) from each 
of the 24 multiple images minus the 9 unknown source 
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positions. Thus each model has ^21 degrees of freedom 
(30 constraints — 9 parameters). The Lenstool model re- 
produces all lensed image positions to an RMS of 1.17'^ 
Assuming a scatter of 1.4^' as explained above, this yields 
a ^ 24 X (l.ir7l.4'0^ = 16.8 with 21 degrees of free- 
dom, for a reduced xl ^ 16.8/21 ^ 0.8. 

The Zitrin spline model with the flat mass profile pre- 
ferring z - 3.5 for MACS0647-JD obtains an RMS of 
1.1^' for a ^ 15, also with 21 degrees of freedom, 
yielding Xu ^ O-'^- When the models are forced to adopt 
2; ^ 11, the best fit is found with a Gaussian-smoothed 
model, yielding an RMS of 2.9", for a ^ 103 and a 
reduced xl ^ 4.9. 

Assuming MACS0647-JD is at z = 11.0, the Lenstool 
model (Fig. 1) estimates magnifications of ~ 8.4, 6.6, 
and 2.8 for JDl, JD2, and JD3, respectively, with uncer- 
tainties of ~ 20%. These magnifications are consistent 
with a F160W = 20 ± 4 nJy (28.2 ± 0.2 mag) source 
magnified by factors of ~ 8.1, 6.8, and 2.1 to match the 
observed F160W fluxes within their ~ 10% uncertainties 
(Table 2). 

LensPerfect models perfectly reproduce all 24 observed 
image positions as input. These data constrain the mass 
distribution and profile well globally but only to a res- 
olution of ~20", or ~130 kpc, roughly the average sep- 
aration between the strongly lensed images. This res- 
olution is insufficient to obtain robust estimates of the 
image magnifications which are strong functions of local 
mass gradients. Here these magnifications are better es- 
timated by adopting priors of light tracing mass as in the 
parametric methods. 

7. PHYSICAL PROPERTIES OF MACS0647-JD 

We estimate the rest-frame UV luminosity and star 
formation rate (§7.1), infer a rough stellar mass (§7.2), 
and place upper limits on the physical size (§7.4). In 
§7.3 we explore SED parameter degeneracies and place 
modest constraints on other properties of MACS0647- 
JD. 

7.1. Rest-frame UV Luminosity; Star Formation Rate 

As described above, we estimate the intrinsic (un- 
lensed) magnitude of MACS0647-JD to be 28.2 ± 0.2 
in F160W. At z ~ 11, the Lyman-a break falls within 
the F160W bandpass, attenuating the observed fiux; the 
rest-frame UV (0.16/im) continuum flux is ~0.25 mag 
brighter. To convert this flux to rest-frame UV absolute 
magnitude Muv,ab^ we add three terms. Most signifl- 
cantly, the magnitude is brighter by the distance mod- 
ulus ~50.3. The flux per unit frequency is also dim- 
mer by a factor of 1 -\- z {^2.7 mag) simply because the 
rest frame samples a higher frequency. We also derive 
a small color term of ^0.1 mag as we switch from the 
blueshifted F160W fllter (-0.13/im) to a tophat fllter 
centered on 0.16//m for comparison with previous mea- 
surements. Combining these terms, we flnd Muy ~ 
— 19.5. Converting this to UV luminosity at a distance 
of 10 pc, we flnd Luv - 2.8 x 10^^ erg s"^ Rz'^. 

This rest-frame UV luminosity can be generated by 
a star formation rate (SFR) of ^AMq yr~^ assuming a 
Salpeter (1955) IMF with mass limits 0.1-100 Mq (Ken- 
nicutt 1998). The ionizing efflciency could be increased 
by a factor of -1.8 for a Chabrier (2003) IMF or -3 for a 
top-heavy IMF (Bruzual & Chariot 2003; Schaerer 2003; 



Stiavelli et al. 2004), with the latter generally realized 
in simulations of high-redshift galaxies (e.g., Abel et al. 
2002; Bromm et al. 2002). Stellar rotations may also in- 
crease this efflciency by a factor of —2-5 (Levesque et al. 
2012). Given these and other uncertainties, a star for- 
mation rate of —1 Mq yr~^ or lower could generate the 
Luv derived for MACS0647-JD. 

This luminosity Luv is or perhaps a few times 
brighter than L*, depending on which extrapolation we 
assume to estimate this characteristic luminosity at z ~ 
11 (Bouwens et al. 2008; Robertson et al. 2010; Bouwens 
et al. 2011a; Bradley et al. 2012b). Based on the esti- 
mated luminosity function (§8) and our lens magniflca- 
tion model, we flnd that a z ~ 11 galaxy lensed to 26th 
magnitude does in fact have an ~80% likelihood of being 
intrinsically brighter than L*. 

7.2. Stellar Mass 

Meaningful observational constraints on the stellar 
mass of MACS0647-JD would require rest-frame optical 
photometry redward of 0.4/im (beyond the Balmer and 
4000 A breaks), or 4.8/im observed. However, we may 
infer a stellar mass estimate as follows. 

Speciflc star formation rates (sSFR) of 2-3 Gyr~^ (that 
is, 2-3 Mq formed per year per 10^ Mq total stellar mass) 
are observed on average for galaxies over a remarkably 
broad range of redshifts (2 ^ z < 7; see e.g.. Stark et al. 
2009; Gonzalez et al. 2010; McLure et al. 2011; Bouwens 
et al. 2012b). If this "plateau" continues out to z ~ 11 
and MACS0647-JD has a typical sSFR of -2 or 3 Gyr-\ 
then this combined with our derived SFR would imply 
a stellar mass on the order of ^IO^Mq. The average 
stellar mass of L* galaxies was ^lO^M© at z ~ 7-8 and 
rose to a few times IO^^Mq hy z ^ 2 (Gonzalez et al. 
2010; Labbe et al. 2010; Finkelstein et al. 2010). Based 
on this trend, we may expect the average stellar mass of 
L* galaxies at z ^ 11 to be less than lO^M©. If this holds 
true for MACS0647-JD, it would suggest a higher sSFR, 
more in line with expectations from simulations which 
are in some tension (but perhaps only mild tension) with 
the observed sSFR plateau (e.g., Khochfar & Silk 2011; 
Dave et al. 2011; Weinmann et al. 2011; Behroozi et al. 
2012). 

We conclude the stellar mass of MACS0647- JD is most 
likely on the order of 10^-10^ M©. The lower end of this 
mass range is more compatible with expectations from 
cosmological simulations and galaxy formation models. 
Based on simulations (e.g., Klypin et al. 2011), we may 
expect to flnd a dark matter halo of virial mass ^lO^^M© 
or so within our search volume of a few times 1000 Mpc^. 
This would comfortably host a galaxy of stellar mass 
10^ Mq or so. a stellar mass of lO^M© would be larger 
than expected. 

7.3. Other SED-based Constraints 

In §5.1 we explored a broad range of galaxy properties 
to rule out lower redshift interlopers with a high degree of 
confldence. In this section we quantify the degeneracies 
in those parameters (see also Pirzkal et al. 2012; Pforr 
et al. 2012) along with the modest constraints we obtain 
on them, assuming that MACS0647-JD is indeed at high 
redshift. 

We modeled the observed photometry using the flexible 
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Figure 15. Constraints on redshift, age, star formation history (exponential scale factor r), metallicity (where Zq ^ 0.019), and dust 
extinction (in V-band magnitudes). This is based on fitting the HST+Spitzer photometry integrated over the three images to the flexible 
stellar population synthesis models of Conroy et al. (2009) and Conroy & Gunn (2010) convolved with an exponential star formation 
history oc exp(— t/r) and with Calzetti (2001) dust extinction added (Ry = 4.05). The star formation history was either decaying 
(r/age = 0.3,1,3), rising (r/age = —0.3, —1,-3), constant star formation rate (CSF; r = oo), or occurring in a single burst at "birth" 
(simple stellar population, SSP; r = 0), with equal likelihoods for all eight possibilities. We assumed a flat linear prior for age (though it's 
plotted as log here) up to the age of the universe (as a function of redshift). Within each panel, the marginalized likelihood is plotted as a 
color map (scaled linearly) and confldence contours of 68% and 95% are overplotted as blue and white lines, respectively. 



stellar population synthesis (FSPS) models from Conroy 
et al. (2009) and Conroy & Gunn (2010). As described in 
§5.1, we convolved their simple stellar population (SSP) 
models with a range of star formation histories, including 
a single early burst, exponentially declining, constant, 
and exponentially rising. Though little dust is expected 
at these redshifts (e.g., Bouwens et al. 2012b), we tested 
this assumption explicitly by adding a variable degree of 
Calzetti (2001) dust extinction. We did not add nebular 
emission lines. Only [Oil] (3727A) redshifted to '^4.5/im 
might be a significant concern, though we only modestly 
detect 4.5/im flux in JD2. We assume a flat linear prior 
in age up to the age of the universe at each redshift. 

Our constraints on redshift, age, SFH (r), metallicity, 
and extinction are shown in Fig. 15 based on the inte- 
grated photometry of all three images. We confirm that 
significant dust extinction is unlikely, with a rest-frame 
F-band extinction of Ay < 0.25 mag (95% confidence). 
Constraints on the other parameters are modest, with 
slight preferences for low metallicity and rising or con- 
tinuous SFH in a maximally old galaxy (limited to ^400 
Myr by the age of the universe). 
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Figure 16. Upper limit on MACS0647-JD half-light radius (de- 
convolved and delensed) compared to observed mean galaxy sizes 
from Oesch et al. (2010a) and Mosleh et al. (2012) extrapolated to 
z ^ 11. The intrinsic scatter in galaxy sizes is a factor of ~ 2 as 
measured for well-studied samples at z ~ 4 (Ferguson et al. 2004). 
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Observed F160W AB Magnitude 



Figure 17. Observed number counts for both lensed (squares) 
and field (circles / arrows) galaxies at z ~ 8, 9.6, and 10.8 as 
a function of observed F160W AB magnitude. CLASH discov- 
eries of candidates at z ^ 10.8 (MACS0647-JD; this work) and 
z ^ 9.6 (MACS1149-JD; Zheng et al. 2012) in ~78 square arcmin- 
utes (17 cluster lensing fields) are plotted as the red and yellow 
squares, respectively, in units of counts per unit magnitude, red- 
shift, and square arcminute. The blue square indicates the single 
robust z ^ 7.5 candidate (A1689-zDl; Bradley et al. 2008) identi- 
fied in 11 non-CLASH cluster fields covering 21 square arcminutes 
(Bouwens et al. 2009). Cyan and orange circles represent counts 
in unlensed fields at z ~ 8 and 10.3, respectively from the BoRG 
survey (Bradley et al. 2012b) as weh as CANDELS, ERS, and 
HUDF09 (Oesch et al. 2012a) including UDFj-39546284 (Bouwens 
et al. 2011a). The arrows give both l-a and 90% confidence upper 
limits from Oesch et al. (2012a). The bottom dashed curves are 
expected counts with uncertainties from the Bradley et al. (2012b) 
z ^ 8 luminosity function and extrapolated to z ~ 9.6 and 10.8 
assuming an evolution of dM"" /dz = 0.30, similar to that found in 
e.g., Bouwens et al. (2008, 2012a) (also see Fig. 18). This yields 
F160W m* ~ 27.0, 27.7, and 28.3 at these redshifts, plotted as a 
small star along each curve. We then simulate the lensing of these 
expected counts in our CLASH WFC3/IR survey area using mass 
models of the 17 clusters and excluding the area (~17%) covered 
by foreground objects. These are given as the upper solid curves 
which are consistent with the CLASH discoveries (red and yellow 
squares). 

7.4. Physical Size 

After correcting for the observed PSF, JDl and JD2 
have observed half-light radii < 0.1'^ or delensed < 0.03^' 
(< 0.1 kpc). Based on extrapolations from lower red- 
shifts (Oesch et al. 2010a; Mosleh et al. 2012), we expect 
an average half-light radius of roughly (ri/2) ~ 0.26 kpc 
for a galaxy with a stehar mass --10^ M©. MACS0647-JD 
is likely somewhat less massive (§7.2). Scatter in galaxy 
sizes is large: ~0.3 dex, or a factor of ~2, as found for 
well-studied samples at 3 < z < 5 (Ferguson et al. 2004). 
So our derived ri/2 ^0-1 kpc is on the small side, though 
not beyond expectations (see Fig. 16). Furthermore, we 
may only be detecting a bright star forming knot in a 
larger galaxy. These knots typically have sizes of ~ 0.1 
kpc as observed in high-redshift (5 < z < 8) lensed galax- 
ies (Franx et al. 1997; Bradley et al. 2008, 2012a). 

8. NUMBER COUNT AND STAR FORMATION RATE 
DENSITIES AT Z > 10 

Our discovery of a 26*^ magnitude z ^ 11 candidate in 
17 cluster lensing fields (^78 square arcminutes) agrees 
with rough expectations given observed luminosity func- 
tions (LFs) at lower redshifts extrapolated to higher red- 
shift and propagated through our lens models. 



The LF at z ~ 8 has recently been robustly constrained 
at both faint (Bouwens et al. 2011b) and bright (Bradley 
et al. 2012b) magnitudes (see also Oesch et al. 2012b). 
Based on the combined HST data from the HUDF09, 
ERS, CANDELS, and BoRG pure parallel fields, Bradley 
et al. (2012b) find the z - 8 LF follows a Schechter 
(1976) function with a normalization (/)* = 4.312 ? x 
Mpc~^, characteristic rest-frame UV absolute magnitude 
M^y = -20.26t^;3^, and faint end slope a = -1.98t^;22. 
These data and LF with uncertainties are plotted in blue 
in Fig. 17. 

We then assumed an evolving LF in which M* varies 
linearly with redshift while ^* and a are fixed. Previous 
work (e.g., Bouwens et al. 2008, 2012b) has shown that 
dM"" /dz ~ 0.3 yields good agreement to data at 4 < z < 
8 (see also Fig. 18). So we first assumed that this holds 
out to 2; ~ 11. 

We convolved this evolving LF through our lens models 
for the 17 CLASH clusters studied in this work, account- 
ing for both the brightening of sources and the reduction 
in search area due to the magnifications. Some of these 
models have been published (Zitrin et al. 2011b, 2012b, a; 
Coe et al. 2012; Zheng et al. 2012; and this work) and 
the rest will be detailed in upcoming work. We applied 
masks to the lensed regions, restricting our search area 
to the WFC3/IR observations and excluding regions cov- 
ered by foreground objects (~17% of the total area) ac- 
cording to our SExtractor segmentation maps. Our total 
search area for 17 clusters is ~78 square arcminutes (as 
observed and lensed). 

The resulting expected lensed number counts for z ~ 8, 
9.6, and 10.8 are plotted in Fig. 17. These are consistent 
with CLASH observed '^26th magnitude number counts 
at ^ - 10.8 (this work) as weh as z - 9.6 (MACS1149- 
JDl; Zheng et al. 2012). 

However, the observed z ~ 10 number counts in the 
field are a factor of ^4 lower than expected based on 
an evolving LF such as this, suggesting a sharp drop 
off in star formation density at these redshifts (Bouwens 
et al. 2011a; Oesch et al. 2012a). To test for such a 
dropoff, we allowed for more (or less) rapid evolution in 
M* at 2; > 8, still as extrapolated from the Bradley et al. 
(2012b) LF at z - 8. We found that MACS0647-JD con- 
strains dM*/dz = 0.30toil, and MACS1149-JD1 con- 
strains dW/dz = 0.62to:33, 

yielding a joint constraint 

of dM^'/dz = OAOtoi7' 
We then integrated these LFs (with uncertainties) 
down to 0.05L*_3 {Muv = -17.7) to obtain star for- 
mation rate densities which can be compared directly 
to previously published estimates (Fig. 18). We found 
SFRD = {1.9tli) X 1O"^M0 yr-i Mpc"^ at z - 11 and 
{1.2l:li) X 1O-^M0 yr-i Mpc"^ at z - 9.6. For consis- 
tency with the other measurements derived and compiled 
by Bouwens et al. (2012b), we assumed a star formation 
rate of one solar mass per year produces a UV (0.14/im) 
luminosity of 8 x 10^^ ergs s~^ Hz~^, as from a Salpeter 
(1955) IMF truncated between 0.1 and 125 Mq. 

Our SFRD estimate z ^ 9.6 is consistent with an 
independent estimate of {l.Sttl) x 1O"^M0 yr"^ Mpc"^ 
based on MACS1149-JD1 presented in Zheng et al. 
(2012). Only 12 clusters were searched in that work. 
To account for the larger volume now searched with- 
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Figure 18. CLASH constraints on the cosmic star formation rate density at 2; ~ 9.6 and 10.8 compared to previous estimates from 
< z < 10 as compiled and corrected for dust extinction by Bouwens et al. (2012b). We adopt the z ^ 8 LF from Bradley et al. (2012b) 
and assume M* evolves with linearly redshift. MACS0647-JD constrains dM* /dz = 0.30+o^g, while MACS1149-JD1 (Zheng et al. 2012) 
constrains dM* j dz = 0.62^0 33, yielding a joint constraint of dM* /dz = OAO'^q'^^. We integrate these LFs down to 0.05L*(2; = 3), or 
Mjjv, 1400 ~ —17.7 AB mag, to obtain rest-frame UV luminosity densities (magenta and green points with the gray shaded region giving 
the joint constraint). The yellow point is an alternative estimate of the z ~ 9.6 SFRD from Zheng et al. (2012) rescaled by 12/17 to account 
for the larger search volume behind 17 clusters. An extrapolation from z ~ 8 assuming dM* /dz = 0.30 is shown as the dashed blue line. 
Points along the bottom blue curve are observed while the upper orange curve is corrected for extinction by Bouwens et al. (2012b). We 
assume a star formation rate of one solar mass per year (left axis) produces a UV (0.14/im) luminosity of 8 x 10^'^ ergs s"-*^ Hz"-*^ (right 
axis), as from a Salpeter (1955) IMF truncated between 0.1 and 125 Mq. A Chabrier (2003) IMF would yield lower SFRD by a factor 
of ~1.8. Previous data are from Schiminovich et al. (2005); Oesch et al. (2010b); Reddy &; Steidel (2009); Bouwens et al. (2007, 2011a); 
Oesch et al. (2012a) as described in Bouwens et al. (2012b, see their Fig. 19). 



out additional ~26th magnitude z ^ 10 candidates, 
we rescaled this estimate by a factor of 12/17, yielding 
(l.Stos) X IQ-^M© yr-i Mpc"^. This value is plotted 
in Fig" 18. 

Our joint constraint from both galaxies on SFRD at 
z > 8 is also plotted as the gray shaded region in Fig. 18. 
Given the large uncertainties, we cannot confidently dis- 
criminate between the trend observed at lower redshifts 
{dM'^/dz ~ 0.30) and the sharp drop off suggested by 
the paucity of z ~ 10 galaxies detected in the field. 

The dominant uncertainty in our SFRD measurement 
is the Poisson uncertainty of '^0.7 dex (a factor of ~5) 
given our single detection. Subdominant uncertainties 
include cosmic variance in 17 independent fields (e.g., 
Trenti & Stiavelli 2008); uncertainties in the lens models 
(e.g., Bradac et al. 2009); and incompleteness in our abil- 
ity to detect ~26th magnitude galaxies (25.4 < F160W 
< 26.4) after masking out areas corresponding to fore- 
ground galaxies. Accounting for overlapping regions in 
the source plane (multiple images) would slightly in- 
crease our derived SFRD. 

9. CONCLUSIONS 

We have discovered a candidate for the earliest galaxy 
yet known at z = 10.71q;4 (95% confidence limits) when 

the universe was 427^21 million years old. The galaxy is 
strongly lensed by MACSJ0647.7+7015 producing three 
magnified images, including two observed at ~26th mag- 
nitude AB (-162 and 136 nJy) in HST WFC3/IR F160W 
imaging (~1.4-1.7/im). The intrinsic (delensed) magni- 
tude is ~20 nJy (mag ~ 28.2) based on our lens models 
for a galaxy at z ~ 11. The unattenuated continuum 
is ^0.25 mag brighter (lensed ~0.2 /iJy in the brightest 



image) . 

Spitzer/IRAC upper limits further support the high 
redshift solution. We tested a broad range of lower 
redshift interlopers, including some previously published 
as high-redshift candidates, and showed that none is 
able to reproduce the observed HST+Spitzer photom- 
etry. Galaxies of known types at z < 9.5 are formally 
ruled out at 7.2-<j, with the next most likely alternative 
being an early type and/or dusty galaxy at z ~ 2.5. Our 
Bayesian priors assume a z ~ 2 galaxy is over 80 times 
more likely than a z ~ 11 galaxy, making our z ^ 11 
claim more conservative than if such a prior were ne- 
glected. For MACS0647-JD to be at z < 9.5, it appears 
it would have to belong to a new class of objects not yet 
observed. 

The discoveries of both MACS0647-JD at z - 10.8 
(this work) and MACS1149-JD1 at ^ - 9.6 (Zheng et al. 
2012) in CLASH observations of 17 clusters to date 
are consistent with extrapolations of luminosity func- 
tions observed at lower redshifts (Bouwens et al. 2012a; 
Bradley et al. 2012b), assuming a linear evolution of M* 
with redshift, and as convolved through our lens mod- 
els. If these extrapolations are valid to z > 10, then 
low luminosity galaxies could have reionized the uni- 
verse (Bouwens et al. 2012a; Kuhlen & Faucher-Giguere 
2012). However these extrapolations are in conflict with 
the paucity of z ~ 10 galaxies discovered in unlensed 
fields, suggesting a rapid buildup in star formation den- 
sity between z ^ 10 and 8 (Bouwens et al. 2011a; Oesch 
et al. 2012a). Our data do not allow us to discriminate 
between these two scenarios, given the large uncertainties 
dominated by the Poisson statistics of these two detec- 
tions. 
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MACS0647-JD is likely close to the characteristic lu- 
minosity for a z ~ 11 galaxy (~l-3 L*), producing a 
few Mq yr~^, with an inferred stellar mass of roughly 
-lO^-lO^M©, and a half-light radius of < 100 pc (de- 
convolved and delensed). This is smaller by a factor of a 
few than the average size expected as extrapolated from 
lower redshifts (Oesch et al. 2010a; Mosleh et al. 2012) 
with an intrinsic scatter in sizes of perhaps a factor of ~2 
(Ferguson et al. 2004). The size of < 100 pc is similar 
to the sizes of bright knots observed in lensed galaxies at 
b < z <S (Franx et al. 1997; Bradley et al. 2008, 2012a). 

Thanks to the magnified views afforded us by gravita- 
tional lensing, this galaxy may be studied further with 
existing and future large telescopes, including the James 
Webb Space Telescope (JWST; Gardner et al. 2006; Sti- 
avelli 2009) and extremely large ground-based telescopes 
constructed in the northern hemisphere. Unfortunately 
due to its high declination of +70, it is not accessible to 
southern telescopes such as the Atacama Large Millime- 
ter Array (ALMA). 

This z ^ 11 candidate approaches the redshift limit 
of galaxies detectable by Hubble's WFC3/IR camera. 
Galaxies at z > 12 would drop out completely of the 
F140W filter and to an increasing degree in F160W until 
z ^ 13, when all the light red ward of Lyman- a would be 
redshifted beyond the observable wavelength range. 

We thank Nor Pirzkal for useful discussions regard- 
ing the WFC3/IR backgrounds and contributing the sin- 
gle emission line galaxy scenario. We thank the ref- 
eree for useful comments which helped us improve the 
manuscript. 
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